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SUMMARY 


This report presents the details and results of an extension of the general-purpose finite clement STructural Anal- 
ysis Routines (STARS) to perform a complete linear aeroelastic and acroscrvoelastic analysis. The earlier version 
of the STARS computer program enabled effective finite element modeling as well as static, vibration, buckling, 
and dynamic response analysis of damped and undamped systems, including those with pre-stressed and spinning 
structures. Additions to the STARS program include aeroelastic modeling for flutter and divergence solutions, and 
hybrid control system augmentation for aeroservoelastic analysis. 

Numerical results of the X-29A aircraft pertaining to vibration, flutter-divergence, and open- and closed-loop 
aeroservoelastic controls analysis are compared to ground vibration, wind-tunnel, and flight-test results. The open- 
and closed-loop aeroservoelastic controls analyses are based on a hybrid formulation representing the interaction of 
structural, aerodynamic, and flight control dynamics. 

NOMENCLATURE 
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A 
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AIC 
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B e 
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C d 

t 

t 

C d ,C l ,Cy 
Cl, Cm, C n 
CPM 


c 

D 

DLM 

dB 


plant dynamics matrix in body-fixed coordinates 

aerodynamic influence coefficient matrix 

coefficient matrices of aerodynamic approximation, (J = 0,1 ... ) 

plant dynamics matrix in inertial frame of reference 

solution to aerodynamic approximations A ; from equation (10) 

body axis accelerations 

six-degree-of-freedom, wind-tunnel based aerodynamic simulation data 

aeroelastic influence coefficients 

aeroservoelastic analysis 

body axis perturbation accelerations 

control influence matrix in body-fixed coordinates 

effective stiffness correction matrix defined in equation (A-26) 

control influence matrix in inertial frame of reference 

wing span 

output state matrix in body-fixed coordinates 
elastic damping matrix 
generalized damping matrix 
C + A] 

coefficients of drag, lift, and side force 

coefficients of roll, pitch, and yaw moment 

constant pressure method 

mean aerodynamic chord 

output control matrix in body-fixed coordinates 

doublet lattice method 

decibel 
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M 
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P, Q, R 
Pi, Qi, R\ 
P, 9, r 

Q 

Q 

Q R, Qi 
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degree of freedom 
column of identity matrices 
force matrix 
finite-element model 
feedback controller matrix 
ground vibration survey 
gravity vector 
aeroelastic damping 

open-loop transfer function matrix of plant and all analog elements 

closed-loop transfer function matrix 

hybrid loop gain matrix 

high frequency crossing (see table 6) 

high gain crossing (see table 6) 

identity matrix 

imaginary part 

v/TT 

elastic stiffness matrix 
generalized stiffness matrix 

K + q\o (see equation (15)) 
effective stiffness matrix 
knots equivalent airspeed 
reduced frequency 
discrete set of reduced frequencies 
low frequency crossing (see table 6) 
low gain crossing (see table 6) 
elastic inertia matrix 
generalized mass matrix 

M + q( £) 2 A 2 (see equation (15)) 
forcing function for elastic dynamics 
generalized forcing function 
total roll, pitch, and yaw rate 
reference (trim) roll, pitch, and yaw rate 
perturbation roll, pitch, and yaw rate 
generalized aerodynamic force matrix 
approximation of Q 
real and imaginary parts of Q 
Q without the rigid air loads 
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T a 

T 

U, V, W 
Ui, V,, W l 
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U, V, w 

VW 

W\B,F\B,... 
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x ; 

x, y , * 
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a 
a i 

ck(w) 
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stacked matrix of Q for a discrete set of reduced frequencies 

generalized force 

displacement vector 

modal displacements at sensor location 

dynamic pressure 

velocity and accelerations 

matrix (see equation (18)) 

real part 

reference input vector 
position vector 
see equations (6) and (7) 
reference area of wing 
STructural Analysis Routines 
Laplace variable 
thrust vector 

inertial to body-fixed coordinate transformation matrices for rigid body states, »' = 1, 2, 3 
(see appendix B) 

inertial to body-fixed coordinate transformation matrices for combined rigid, elastic, and 
aerodynamic lag states 

sensor interpolation matrix 

sample time for digital controller 

total velocities in body-fixed frame of reference 

reference (trim) velocities in body-fixed frame of reference 

system input vector 

perturbation velocities in body-fixed frame of reference 
virtual work 

vibration modes (see tables 1 and 2) 
system state vector in body-fixed reference frame 
system state vector in inertial frame 
;th lag state vector 

longitudinal, lateral, and directional displacement 

system output vector 

angle of attack 

trim angle of attack 

equivalent downwash 

angle of sideslip 

aerodynamic lag terms 

inertial to body-fixed coordinate transformation matrices, i = 1, 2 (see appendix B) 
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control surface deflection 

V 

generalized coordinate vector 

0, <t>, 'F 

total pitch, roll, and yaw angles 

01, Oi, 'F, 

reference (trim) pitch, roll, and yaw angles 

o, 4>, 

perturbation pitch, roll, and yaw angles 

A , A , A,', Aj 

inertial to body-fixed coordinate transformation matrices, t = 1,2 (see appendix B) 

r 

system time delay 

Q, Cl 

angular velocity matrices 

Ul 

frequency 

c 

damping 


Subscripts 

B 

body-fixed frame of reference 

D 

drag 

e 

elastic 

I 

inertial frame of reference 

L 

lift 

Q 

pitch rate 

R 

real part of a complex number (inertial frame of reference in appendix B) 

r 

rigid 

x,y,z 

coordinate system 


Superscripts 

T 

matrix transpose 

-1 

matrix inverse 


INTRODUCTION 

Aircraft systems are becoming increasingly complex because of the integration of distinct technologies used 
to attain numerous objectives in the areas of performance, control, flying qualities, maneuver techniques, fuel ef- 
ficiency, and various mission requirements. The design process must respond to specifications from all disciplines 
to achieve these diverse goals and integrate accordingly. Design procedures must account for conflicting objectives 
and the interaction of dynamics from the control system, structure, and aerodynamics. Coupling between these dy- 
namic elements of the model can be treated passively with structural modifications and passive filtering, or actively 
with control mechanisms driven by appropriate control laws. Analysis of the consequences of the design is essential 
to perform safe and effective mission tasks. Aeroservoelastic analysis is used to investigate the potential problems 
arising from high bandwidth control of relatively flexible aircraft by combining linear models of structure, unsteady 
aerodynamics, and control system into one dynamic system. The models are augmented to address stability and 

performance issues. 
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Vibration, aeroelastic, and aeroservoelastic analyses are presented in this report, computed with the STructural 
Analysis Routines (STARS) program, an integrated analysis tool evolved from an earlier (ref. 1) finite element 
structural analysis program. A number of available computer programs perform only individual aspects of these 
analyses (refs. 2-5). The STARS program was created to analyze complex practical systems such as aircraft and 
is designed as a unified, compact, modular, and highly graphics oriented analysis tool. Extensive verification and 
application of the program has been made at the NASA Dryden Flight Research Facility in support of flight programs, 
such as the X-29A, through correlation of analysis with ground-test and flight-test data. The formulation of the 
mathematical procedures, a program outline and description, and selected results of these analyses are presented. 

FORMULATION OF NUMERICAL PROCEDURE 

Extensions made to the STARS program include aeroelastic modeling and control system augmentation to the 
aero-structural dynamics. An approximation for the unsteady aerodynamics is performed with a Padd fit to the 
aerodynamic influence coefficients (AIC) to generate a state-space aero-structural dynamic model. This model is 
transformed from the inertial to body-axis coordinate system for control system augmentation. Aeroservoelastic 
analysis is achieved for either analog or digital controllers with hybrid frequency responses and eigenvalue solutions 
for closed-loop modal behavior. 

The matrix equation of motion for structures relevant to the current analysis has the following form 

Mq + Cq + Kq + q\ e (k)q = P(t) (1) 


in which 

M is the inertia matrix, 

Q the damping matrix, 

K the elastic stiffness matrix, 

q dynamic pressure, 

k reduced frequency jy, w and c being the natural frequency and mean chord length, respectively, 

A e ( k) the aerodynamic influence coefficient matrix calculated for a given 

Mach number and set of k values, 
q the displacement vector, and 

P (t) the external forcing function. 

The free vibration solution is first affected (refs. 1 and 6) on the matrix formulation 

Mq + Kq = 0 (2) 

yielding the frequencies, w, and mode shapes <1>. Applying a transformation 

q = <J>i? (3) 


to equation (1) and pre-multiplying both sides by O r yields the generalized equation of motion 

Mrj+ trj+ ftr? + gQ(fc)r?= P(t) (4) 

where M = ft = Q = <b T 'A e <I>, P = the modal matrix <I> — [<t> r <I> e <I> 5 ] , 

and the generalized coordinate tj = [7? r rj e TjJ , thereby incorporating rigid body, elastic, and control surface motions, 
respectively. A general Lagrangian formulation of the equations of motion for an aircraft maneuver analysis is given 
in appendix A. 
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The generalized aerodynamic force matrix Q(fc) may be approximated (ref. 7) with Padd polynomials in 
i*k(= i*fy = jy, where the Laplace variable s = i*w) 


0 (A:) = A 0 + + (i*fc ) 2 A 2 + rrr — — A 3 + — — — A 4 + 


i*A: + /?i 


i*k + P 2 


(5) 


with aerodynamic lag terms P } (assume j = 1 , 2 ), and 

i*k k z 


_ ^ i'kpj 

\*k + (3j k 2 + 0j ^ k 2 + P 2 

The rigid airload coefficients assume the following form 


Ao = Qr( k \ ) 

Q/(*l) ^3_ A^ 

A:i Pi P2 


A, = 


where k\ is the smallest reduced frequency, with a value near zero, used to compute A ; for j — 0, 1 , 2, . . . Separating 
real and imaginary parts in equation (5) yields 


Qfl(fc) 1 = ■Qdk ) )- -Ao 

= [^ 2 ' pV P75T 1 ] 
= s R (k)\ 

Q j(fc) 4 


A 2 

a 3 

A4 


( 6 ) 


Qi(fc) = 


- [ 0 T&ff 1 ] 

= S/(Jfc)A 


A 2 

A 3 

A 4 


(7) 


The unknown coefficients A 3 and A 4 can be determined by substituting the previous expression for A] into 
equation (7). However, the resulting solution is sensitive to the choice of pj for approximating rigid airloads. If the 
elements of the aerodynamic damping matrix Aj are replaced with known damping coefficients (steady aerodynamic 
derivatives), then the solution for rigid airloads becomes insensitive to the Pj values. 

For a chosen number of values of reduced frequencies fc,, equations ( 6 ) and (7) may be combined as 


' Q R (k 2 ) 


' S R (k 2 ) 

Oi(fe) 


S/(fc 2 > 

Qr( kuF) 


Sh( kfiTF- 1 ) 

( h(kNF ) 


r 

c n 

/—v 

1 


A 2 

a 3 

A 4 


( 8 ) 


or 


0 = SA 


(9) 
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and a least-square solution 


A = [S r S]-'S r 0 (10) 

yields the required coefficients A2 , A3 , A4 . The procedure is easily extended for a larger number of lag terms, if 
desired. Equation (4) may then be rewritten, assuming simple harmonic motion, as 



A A A sc { sc \ ^ 

M77 + Cr]+ Ktj + g[A 0 rj + A] —77 + A 2 ( — J t) + A 3 Xi + A 4 x 2 + ■••] = 0 

( 11 ) 

such that 

ST) 

Xj ~ s+ 2 -fPj 

(12) 

from which 

(2V\ „ 



X; + (^ — ) PjXj = 17 

(13) 

Collecting like terms gives 



(K + q\o)rj+ (C + g^yA 1) 77+ MVI + q A 2 ^ rj+ qAj\i + gA 4 x 2 + •• 


= 0 


or 


Krj+C77+Mrj + gA3 Xj + gA 4 x 2 + ■ • - = 0 
Rewriting equations (13) and (15) as one matrix equation 

I 


or 


and 


M 





0 

1 


V 


-k 

-t 


Xl 


0 

1 


_ *2 . 


0 

1 


0 


-t&I -I 



V 


V 


Xl 


_ X 2 


M'x'= K'x' 

x' = (M') -1 K'x' 
= Rx' 


Now rearranging the state-space vector x' as 


(14) 

(15) 


(16) 


(17) 


(18) 


x " = [(VrVA'O'XiXz) (%%)} 

= L* U J 


(19) 


equation (18) may be partitioned as 


X 


Ri.i 

Ri ,2 


X 

u 


1*2,1 

1*2,2 


u 


( 20 ) 


where the first matrix equation denotes the plant dynamics and the second represents the dynamics of control modes. 
In the case of plant dynamics, the state-space equations become 


x = Ax + Bu 


( 21 ) 
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the associate matrices and vectors being defined as 


A plant dynamics matrix 

B control surface influence matrix 

x generalized coordinates in inertial frame 

u control surface motion input into plant 

and in which the terms Ax and Bu represent the airplane dynamics and forcing function due to control surface 
motion, respectively. 

Coordinate Transformation 

To incorporate control laws designed to control body-axis motions, it is necessary to transform equation (21) 
from the Earth-fixed (inertial) to the body-fixed coordinate system. Since no transformations are applied to elastic 
and aerodynamic lag state vectors, a transformation of the form 

x = — T 3 )x + t 2 -'6u 

= Ax + Bu (22) 

in which 

ti = 

Coordinate transformation matrix Ti yields the required state-space equation in the body coordinate system. A 
detailed description of the transformation procedure is given in appendix B. 

Determination of Sensor Outputs 

The structural nodal displacements are related to the generalized coordinates by equation (3) and the related 
sensor motion can be expressed as 

q 4 = T 4 «t>T> 

= C 0 x (23) 

where C 0 = [T 4 <h000], and in which J s is an interpolation matrix. Similar relations may be expressed for sensor 
velocities and acceleration as 

Qa 1 _ T T 3 Oi7 

Qs J [ 

(24) 

= Cjx (25) 

where 

T a «h 0 0 0 

Cl = 

0 T a O 0 0 

Pre-multiplying equation (22) by Ci 

Cix = CiAx + CiBu 

= C 2 X + D 2 U (26) 
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Adjoining equations (23) and (25) 


or 



Co 

C 2 


x + 



(27) 


y = Cx + Du 


(28) 


which is the required sensor output relationship signifying motion at the sensors due to body motion (C ) and control 
surface motions (D ). 


Augmentation of Analog Elements and Controller 


Equations (22) and (27) represent the complete state-space formulation for the aircraft incorporating structural 
and aeroelastic effects. To conduct an aeroscrvoelastic analysis, it is essential to augment such a formulation with 
associated analog elements such as actuators, sensors, notches, and pre-filters along with the controller. Denoting 
the state-space equation of one such typical element in series as follows 

x (0 = A ( *V 0 + B (,) u (0 

y<*> = C ( ‘>x (<) + D (0 u (<) (29) 


which can be augmented to the plant equations (22) and (27) as 


X 


*(••) 



A 0 

A <9 


x 

x<‘> 


B 

B<’>D 


[u] 


or 


*(0 = A(,)X ( i) + B(,)U 


noting that u (,) = y from equation (27). Also 


y 


C 0 

X 


D 

y<* 


D (, '>C C ( ‘) 

x (0 

+ 

D (l >D 


[u] 


or 


y<0 = C(i)X(i) + D (1 )U 


becomes the new sensor output expression. 


(30) 

(31) 


(32) 


All analog elements, including a controller, can be augmented similarly at the input and the output of the plant. 
Figure 1 shows a typical feedback control system with controller G . For such a system the three sets of relevant 
matrix equations are 


x = Ax + Bu 
y = Cx + Du 

u = r - Gy (33) 

where equation (32) is the feedback equation. The required transfer functions may be obtained by Laplace transformation 

sx(s) = Ax(s) + Bu(a) 
y(s) = Cx(s) + Du(s) 

u(s) — r(s) — G(s)y(s) (34) 
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From equation (33) 


x(s) = [al -A] _, Bu(a) (35) 

and substituting equation (34) in equation (33) yields the required open-loop frequency response for zero 
initial conditions 


y(s) = [C(«I — A)-'B + D]u(s) 
= H(a)u(s) 


(36) 


where H (a) is the open-loop transfer function without the controller. To obtain the closed-loop frequency response, 
equation (35) is first substituted in equation (33) to get 


u(s) = r(s) — G(s)H(a)u(s) 

= [I + G(s)H(a)] _1 r(s) (37) 


and then using equation (36) gives 

y(s) = H(s)[I + G(s)H(s)] -, r(s) 

= H(a)r(s) (38) 

in which H ( s) is the closed-loop transfer function. The frequency response plots can be obtained from the transfer 
matrices H ( s) or H ( s) as the case may be. Associated damping and frequency values of the system may also be 
calculated by solving the eigenvalue problem of the relevant A matrix for various fc, values or dynamic pressures 
and observing the changes in sign of the real part of an eigenvalue to detect instabilities. 

In the case of a digital controller, a hybrid approach is adopted for the frequency response solution. Thus if A 
B\ C', and D' are the state-space matrices associated with the controller, the related transfer function is simply 
given by 

G(z) = C'[zl — A'r'B' + D' (39) 

and the frequency response formulation (ref. 8) for the hybrid analog-digital system with time delay r and sample 
time T can be written as 

, , , fU(s)[ZOH]\ , , 

y (*) = G(z)[, =e . r] f Ju(s) 

= G(s,T)H*(s,t,T)u(s) ( 40 ) 


in which 

H(s) is the plant transfer function with all analog elements 

[ZOH] is the zero-order hold complex expression (= 

and where H*(s,r, T) is now the loop gain transfer function of the hybrid system. The closed-loop frequency 
response may be obtained as before by using equations (33) and (39) 

y(s) = |h *( s, t , T") [ I + G(s,T)H*(s,T,T)r 1 } r(a) 

= n*(s,T,T)r(a) (41) 

To calculate the damping and frequencies, modes with natural frequencies much beyond the Nyquist frequency 
are truncated. The analog plant dynamics matrix A is then transformed into the z-plane by standard discretization 
procedures and augmented to controller dynamics A '. Appropriate eigenproblem solution of the final matrix yields 
the required results, as previously discussed. 
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DESCRIPTION OF THE STARS PROGRAM 


The STARS program primary modules are shown in figure 2. The associated schematic for the aeroservoelastic 
analysis is shown in figure 3. Some highlights of each module are presented herein. 


Preprocessor 

This module enables automatic generation of finite element mesh and associated STARS input data for any con- 
tinuum. The interactive graphics program is capable of generating complex structural systems through interpolation, 
up ication, mirror-imaging, and cross-sectioning of modular representative forms using menu or command options. 
This program is designed to run on a wide variety of commonly used graphics terminals. 


Finite Element Modeling 

This module is made up of many finite elements which are suitable for modeling complex practical structures. 
Some basic elements are as follows 


1. One-dimensional rods, bars, and beam elements, 

2. two-dimensional, triangular and quadrilateral elements pertaining to membranes, shear, plate bending, and 
shells including sandwich and composites, 

3. three-dimensional, tetrahedron, prism, and hexahedron elements, and 

4. one- and two-dimensional finite dynamic elements. 


Die desirable features of this module include random data input, bandwidth minimization, multiple coordinate 
systems, mechanical and thermal loading, general interdependent deflection boundary conditions, automatic node, 
and element generation as well as plot of geometry. 


Aeroelastic Analysis 

for Ih neWly devel0 fl ped COnstant pressure method < CPM > (ref- 9) is used for computing unsteady aerodynamic forces 
for the supersonic flow, whereas the doublet lattice method (DLM) (ref. 5) is used for the corresponding subsonic 
case. Both k and p-k aerodynamic stability (flutter and divergence) solution methods are available in this module. 

Aeroservoelastic (ASE) Controls Analysis 

In this module, the aerostructural problem is recast in the Laplace domain and the unsteady generalized aerody- 
namic forces are curve fitted, using Pade and least-square approximations, to compute state-space matrices. Such 
matrices are augmented with analog elements such as actuators, sensors, prefilters, and notch filters. 

F ° r a " ana, °g controller, the state-space matrices are appropriately augmented into that of the entire system and 
the related transfer function produces the open- and closed-loop frequency responses as desired. The associated 

modal frequency and damping values may also be derived by solving the eigenvalue problem of the relevant state- 
space (A) matrix involving the plant dynamics. 

In the case of a digital controller, a hybrid approach is adopted in which the associated transfer function is suitably 
combined with that of the original analog system. This combination yields the hybrid transfer function from which 
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the open- and closed-loop frequency responses may be computed as desired. To obtain the modal frequency and 
damping values, the analog A matrix is first discretized and added to the corresponding matrix for the controller. 

Furthermore, stability analyses for open-loop flutter and divergence or closed-loop control augmentation aero- 
servoelastic analysis (ASE) can be performed by solving the appropriate eigenvalue problem. The former analysis 
can be compared with the aeroelastic analysis using k and p-k methods, whereas the latter analysis proves to be 
useful for comparing relevant flight-test results of such unstable aircraft as the X-29A. 


Postprocessor 


Extensive interactive plotting facilities are available for solution results pertaining to each module of the STARS 
program. These facilities include contour lines of deformations and stresses as well as time-dependent funcUons of 
dynamic responses, mode shapes, flutter and divergence plots for k, p-k and root locus methods, and p ase g 
plots pertaining to the ASE analyses as well as the corresponding damping and frequency plots. 

NUMERICAL RESULTS 

The presently developed integrated aero-structural-control analysis program STARS has been used extensively 
to solve related problems of the X-29A forward-swept wing research aircraft (refs. 10 and 1 1) (fig. (4)) The X-29A 
aircraft has thin wings of composite materials that are aeroelastically tailored to eliminate structural divergence 
within the flight envelope. Full-span, double hinged, variable-camber flaperons and shake flaps operate with fu 
authority variable-incidence canards to yield minimum trim drag. The presence of the supercritical airfoil enables 
efficient transonic cruise performance and high transonic maneuvering. The canard configured aircraft is upto35- 
percent statically unstable, thus requiring appropriate feedback controls for augmented static stability, 
these combined technologies result in significant improvement in overall aerodynamic and structural performa c , 
they may also cause adverse dynamic interaction of the flight controls with the flexible structure if not integrated 
properly. Therefore an ASE analysis assumes a very important role in the design process. 

In this report a comparison of ASE as well as vibration and aeroelastic analysis of the X-29A aircraft between 
the STARS calculated results and test results is presented in detail. 


Structural Dynamics and Aeroelastic Analysis 

A finite element model of the X-29A aircraft is shown in figure 5. This somewhat reduced dynamic model 
(aooroximately 3000 dof) was derived from the detailed stress model (7100 dof) by using an equivalent shell method 
for the various appendages to reduce free vibration analysis time. The fuselage was modeled with line elements and 
the vibration problem for the entire structure was solved by a block Lanczos procedure (refs. 1 and 12). A summa^ 
of such free vibration analysis results is compared with those obtained by the ground vibration survey (CVS as 
shown in tables 1 and 2 for the symmetric and antisymmetric cases, respectively. The corresponding node line 
comparisons are shown in figures 6 and 7. A more detailed presentation of the vibration results can be found 

reference 13. 

Extensive aeroelastic analyses were performed to determine flutter and divergence speeds of the aircraft em- 
ploying the k, p-k, and root-contour methods for comparison purposes; associated STARS unsteady aeK>dynanuc 
paneling is shown in figure 8. Such results from the STARS model are compared with GVS results in tables 3 and 4, 
where the column labeled STARS (ASE) represents the results of the eigensolution of the plant A matrix, whereas 
figures 9 through 12 show the various flutter solution plots. 

In general, good correlation of results for various cases was observed from the results presented herein. A more 
detailed mode-by-mode comparison of these results is presented in reference 13. 
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AEROSERVOELASTIC ANALYSIS-COMPARISONS WITH MEASURED DATA 


Hie [ A , B , C , D] dynamic system computed by STARS is compared with measured data in the form of a lin- 
earized aerodynamic database for the X-29A aircraft, called AER07, and flight data. The AER07 is a nonlinear, 
six dof, rigid body wind-tunnel database with flexibility corrections to the rigid data, heretofore denoted as flexibi- 
lized data. The comparisons consist of 

1. Open-loop roots: the rigid body modal characteristics from linearized AER07 data, the STARS rigid body 
modes, and the STARS model with all rigid and elastic modes included, are represented. 

2. Frequency responses: loop gain and longitudinal closed-loop responses are computed with measured GVS 
and STARS modal data aircraft dynamics augmented to a digital flight control system. Stability margins are 
also compared. 

3. Flight data responses: longitudinal loop gains computed from STARS are compared with those generated from 
flight data and linearized AER07 data, and 

4. Flight measured modal data: flight determined modal damping and frequency results are compared to STARS 
results for the first two symmetric modes and first three antisymmetric modes. 

The four flight conditions analyzed are 

1. Mach 0.7 at 40,000-ft altitude, 

2. Mach 0.9 at 30,000-ft altitude, 

3. Mach 0.9 at 15,000-ft altitude, and 

4. Mach 1.2 at 30,000-ft altitude. 

Condition 1 is a low dynamic pressure case, 2 is the X-29A aircraft design point, 3 is a high dynamic pressure 
condition, and 4 is a supersonic point where the CPM was used to calculate the theoretical data. 


Open-Loop Roots 

The rigid body roots are tabulated in tables 5(a), (b), (c), and (d) for the longitudinal modes (phugoid and short 
period) and lateral-directional modes (spiral, roll, and dutch roll). Complex poles are written as (w, Q. The STARS 
roots are compared to linearized rigid (without flexibility corrections) and flexibilized AER07 data. Tne STARS 
rigid body model (rigid) contains all unsteady lag states on the rigid body and control modes, but no elastic modes or 
AER07 flexibilized data corrections to the generalized aerodynamic stiffness and damping matrices (appendix C) 
are included. The models generated with lag states on only rigid and control modes by incorporating AER07 flexibi- 
lized stability derivatives into the generalized aerodynamic stiffness and damping matrices are designated as Rigid- 
flexibilized in the tables. The flexible model (flexible) consists of all rigid, elastic, and control modes with lag states 
included. Flexible models are corrected with rigid AER07 derivatives (not flexibilized) since the flexibility effects 
are already modeled by incorporation of the elastic modes. 

The columns labeled Full contain roots computed from a full-order model, either containing only the rigid modes 
(rigid and rigid/flexibilized) or both rigid and elastic modes (flexible). Columns designated Reduced represent roots 
computed from a reduced form of the full-order model by only including the static contribution of the excluded 
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states Those states describing rigid body dynamics are retained and all other states are excluded in the reduced- 
order model to represent the static contribution of flexibility as well as unsteady aerodynamics on the rigid model. 
Hence, the state vector is partitioned into rigid body states (r), and all other states (o) in the model, as 



and the A , B , C system is similarly partitioned as 



where D is unchanged. Then the reduced model is given by 

Xr = (A rr - AroA^A^X,. 

+ ^B r — AqtAqq B d ) u 
y — (C r — CoAj Aqt^ x r 
+ (D - C 0 A“'B 0 )u 

Residualization of the states in this way results in a reduced-order model reflecting the effects of the deleted elastic 
as well as unsteady dynamics on the rigid body modes. The coupling may be significant especially when the aircraft 

is relatively flexible (ref. 14). 

The root comparisons between the STARS corrected models and AER07 models are quite good for all modes 
except the phugoid and spiral modes. In fact, STARS predicts totally damped phugoid characteristics instead of the 
expected oscillatory behavior. These very low frequency modes are difficult to predict with a Pad6 fit over a wide 
reduced frequency range. Results using the generalized forces at one or two of the lowest reduced frequencies with 
zero or first-order Pad6 fits show more reasonable behavior for both phugoid and spiral modes. This procedure is 
recommended if aeroelastic effects are being investigated on only these low frequency modes. The other modes with 
more reduced frequencies included, however, are reproduced very well with the present STARS implementation of 

AER07 data. 

The uncorrected STARS model does not have a well defined phugoid because the forces in the forward velocity 
dof are zero (appendix C). Hence, a root at zero is evident in all cases. The other root is therefore meaningless. 
A comparison of the rigid AER07 model and the uncorrected rigid STARS model shows how well the theoretical 
STARS aerodynamic model compares to a wind-tunnel derived model. The roll and dutch roll modes demonstrate 
excellent agreement, while the short period is not as close in agreement, but reasonable. 

Aeroelastic effects are most noticeable at flight condition 3, the high dynamic pressure transonic condition, where 
differences between rigid AER07 and flexibilized AER07 are most evident, but still not significant. The STARS 
program also demonstrates flexible effects when comparing the rigid roots with flexible model roots. Similar minor 
differences between full-order and reduced-order solutions are also evident, showing that flexibility effects on rigid 
body stability are essentially negligible. 

In general, the comparisons validate the ability of STARS to reproduce rigid body dynamics with 
flexibility effects. 

Frequency Responses 

Loop gains and longitudinal closed-loop frequency responses were computed and stability margins determined 
for the four flight conditions with the digital normal mode (ref. 11) control system augmented to the aircraft states. 
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The longitudinal and lateral-directional control system block diagrams are shown in figures 13 and 14, respectively. 
Hybrid frequency responses (ref. 8) were calculated to account for aliasing effects of dynamics with frequencies 
higher than the Nyquist frequency of 20 Hz (126 rad/sec). Points designated OL in the diagrams represent open- 
loop response locations. The Cl locations in figure 13 are the longitudinal closed-loop input-output signals from 
filtered pilot stick position to aircraft pitch rate response since pitch rate is the dominant feedback variable. 

Stability margins for all augmented models are displayed in tables 6(a), (b), (c), and (d). Loop gains are plotted in 
figures 15 to 18 (longitudinal), figures 19 to 22 (lateral), and figures 23 to 26 (directional) for all four flight conditions 
analyzed, using measured GVS modal data and modes generated from the STARS. Aircraft dynamics include all 
rigid, elastic, and control modes with lag states and AER07 flexibilizcd data incorporated. Control system dynamics 
include all elements shown in figures 13 and 14. 

Longitudinal frequency response overplots in figures 15 to 18 show some dissimilarity between 30 and 100 rad/sec 
due to the first two symmetric elastic modes. Pitch rate response from the first wing bending and first fuselage ver- 
tical bending modes is more pronounced from the measured GVS modal data than from STARS modes. Notch 
filters were designed for the X-29A aircraft to gain stabilize all modal responses to less than -6 dB in the loop gains 
throughout the flight envelope. Condition 3 at Mach 0.9, 30,000-ft altitude was one of the most critical points in 
the longitudinal axis. Lateral responses in figures 19 to 22 show fairly benign modal contributions to feedbacks, 
and the measured model is not dramatically different from the STARS model. Supersonic condition 4 shows some 
offset between the gain curves due to differences in aileron effectiveness. The directional responses in figures 23 to 
26 demonstrate a much more significant discrepancy between 50 and 100 rad/sec, amounting to 20 to 30 dB in this 
frequency range. Modal contributions from the first three antisymmetric elastic modes to the yaw rate feedback are 
10 to 30 times greater from the measured GVS data than from the STARS model. 

These observations demonstrate an important aspect of the ASE modeling process. A critical issue for an accurate 
analysis is the sensor feedback motion caused by modal dynamics. Since the feedback sensors are usually located 
in the fuselage, the fuselage motion at the sensor locations must be determined precisely. Sufficient accuracy is 
difficult to obtain from a finite element model of the fuselage. A very detailed high-order stress model is needed to 
obtain the proper sensor motion due to fuselage dynamics. Slight differences in mode shapes between the model and 
actual aircraft can result in significant discrepancies in feedback dynamics between the model prediction and actual 
aircraft. A ground test is essential to determine the actual feedback due to modal dynamics. The X-29A aircraft 
ground tests demonstrated that the discrepancies were extremely significant (ref. 4). 

Longitudinal closed-loop plots in figures 27 to 30 for all four flight conditions reveal the flexibility of the GVS 
aircraft model compared to the STARS model. Closed-loop models are significantly different representations of the 
dynamic interaction between the digital control system and flexible-unsteady effects, in contrast to the equivalent 
open-loop description, since the dynamics beyond the Nyquist frequency are propagated through all discrete and 
continuous elements simultaneously. The consequences of sampling and aliasing on these dynamics contribute to 
the modal interaction between the controller and aircraft model. Analysis using the entire closed-loop spectrum 
may expose possible problems in the feedback mechanism due to the interface of discrete and continuous dynamics. 
The relatively small differences between STARS and GVS models in the longitudinal loop gains of figures 15 to 18 
are magnified when the loops are closed in the corresponding closed-loop plots of figures 27 to 30, yet closed-loop 
modal damping is adequate as a result of notch filter design. In these cases, closed-loop plots seem to exaggerate 
minor differences in modal damping between GVS and STARS models, which also depend on phase relationships 
between modal dynamics and feedback sensor motions. 

Flight Data Responses 

Loop gains were computed for the longitudinal axis from flight data (ref. 15) and are compared with the aug- 
mented flexible (including all rigid, elastic, control, and lag states) STARS models in figures 31 to 34 for all flight 
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conditions analyzed. Correlation is excellent for the spectrum range of the measured data. The AER07 models are 
overplotted for reference. The STARS models match the flight data as well as AER07 models for all flight conditions 
shown. These STARS models are corrected with rigid AER07 coefficients as outlined in appendix C. Hence, the 
plots show that STARS models fiexibilized by incorporating the calculated mode shapes properly represent the rigid 
body motion. 


Flight Measured Modal Data 


Figures 35 and 36 portray good agreement between flight measured modal results and STARS calculations for the 
symmetric wing first bending and fuselage first vertical bending modes. The STARS data are obtained by calculating 
the roots of the flexible model augmented to the digital control system. This procedure is possible only for those 
modes with natural frequencies up to and near the Nyquist frequency. Such similarities between STARS and flight 
results validate the STARS formulation for the entire flight envelope. Antisymmetric correlations shown in figures 37 
to 39 for the antisymmetric wing first bending, fuselage first lateral bending, and fin first bending also show the 
effectiveness of STARS to predict elastic behavior. The only noticeable difference between the STARS results 
and flight data is in the antisymmetric wing first bending damping and frequency at lower altitudes. The STARS 
program predicts a higher damping and lower frequency and is therefore less conservative at these conditions. At 
higher altitudes the correlation is much better for this mode. 


CONCLUDING REMARKS 


The mathematical formulation of the procedures to perform vibration, aeroelastic, and aeroservoelastic analy- 
ses with the structural analysis routines program has been presented, along with numerical results for the X-29A 
aircraft. Interaction of structural dynamics, aerodynamics, and flight control systems is modeled using both a finite 
element model and measured modes to predict flutter-divergence and aeroservoelastic phenomena. Correlation with 
measured data from ground tests as well as flight measured results is used to validate the capability of STARS to 
analyze complex dynamical systems such as flexible aircraft controlled with digital control systems. 


Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, March 20, 1989. 
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Table 1. Free vibration analysis results for X-29A symmetric half-aircraft. 


Primary 

motion 

STARS 

frequency, 

Hz 

GVS 

frequency, 

Hz 

Rigid body 

0.00 

n/a 

Wing first bending (W 1 B) 

8.96 

8.61 

Fuselage first vertical bending (FIB) 

12.87 

11.65 

Fuselage second vertical bending ( F2B ) 

19.03 

24.30 

Canard pitch (CP) 

21.02 

21.70 

Wing second bending ( W2B ) 

26.28 

26.30 

Wing first torsion (WIT) 

30.30 

36.70 

Canard bending pitch (CBP) 

47.70 

42.20 

Wing third bending (W3B) 

49.52 

51.50 


Table 2. Free vibration analysis results for X-29A antisymmetric half-aircraft. 


Primary 

Motion 

STARS 

frequency, 

Hz 

GVS 

frequency, 

Hz 

Rigid body 

0.00 

n/a 

Wing first bending (W\B) 

10.08 

11.3 

Fuselage first bending (FIB) 

12.35 

12.5 

Fin first bending (Fin IB) 

17.18 

15.2 

Canard pitch (CP) 

21.52 

21.9 

Wing first torsion (W IT) 

27.15 

26.8 

Wing second bending (W2B) 

32.88 

34.8 

Fin second bending (Fin2B) 

41.58 

45.2 

Wing third bending (W3B) 

45.85 

51.7 

Fin first torsion (Fin IT 1 ) 

48.95 

50.0 

Inboard flap torsion (IFT) 

50.83 

51.0 


Table 3. X-29A aircraft flutter and divergence speeds for symmetric modes. 


Mode 

Divergence point, KEAS 

Flutter point, KEAS 

STARS 

k 

STARS 

ASE 

GVS 

STARS 

k 

STARS 

ASE 

GVS 

W\B 

833 

819 

865 

— 

— 



CP 

912 

910 

1017 







CBP 

— 

— 

— 

— 

— 

694 

W3B 

— 

— 

--- 

1275 

1231 

1222 
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Table 4. X-29A aircraft flutter and divergence speeds for antisymmetric modes. 


Divergence point, KEAS Flutter point, KEAS 


Mode 

STARS 

k 

STARS 

ASE 

GVS 

STARS 

k 

STARS 

ASE 

GVS 

W\B 

F\B 

838 

833 

808 

848 

797 

924 

CP 

W2B 

913 

918 

980 

1143 

1157 

1315 


Table 5. Comparison of rigid body modal characteristics between STARS 
and linearized aerodynamic database for X-29A aircraft. 


Model 


Mode 

Linear 

rigid 

AER07 

STARS rigid 

Linear 

flexibilized 

STARS rigid/flexibilized 

STARS flexible 

Full 

Reduced 

~ AER07 

Full 

Reduced 

Full 

Reduced 

(a) Mach = 0.7, altitude = 40,000 ft. 

Phugoid 

(0.7,0.08) 

0.0,0.11 

0.0,0.12 

(0.07,0.08) 

-0.01,-0.12 

-0.001,-0.12 

-0.001,-0.12 

-0.001,-0.12 

Short period 

2.1, -2.7 

0.9 -1.3 

0.9,- 1.3 

2.1,— 2.7 

2.3,— 2.7 

2.3, -2.6 

2.3, -2.7 

2. 3, -2.6 

Spiral 

0.03 

0.0005 

0.0005 

0.03 

0.5 

0.003 

0.5 

-0.001 

Roll 

-1.1 

-1.1 

-1.3 

-1.1 

-1.3 

-0.9 

-1.3 

-1.0 

Dutch roll 

(1.9,0.05) 

(1.6,0.07) 

(1.6,0.07) 

(1.9,0.05) 

(2.2,0.05) 

(2.0,0.08) 

(2.2,0.05) 

(2.0,0.08) 

(b) Mach = 0.9, altitude = 30,000 ft 

Phugoid 

(0.08,0.17) 

0.0,0.32 

0.0,0.35 

(0.08,0.17) 

-0.003,-0.28 

-0.003,-0.28 

-0.003,-0.31 

-0.003,-0.32 

Short period 

3.2,— 4.5 

1.1, -2.3 

1.1, -2.3 

3.3, -4.6 

3.6, -4.6 

3.5, -4.4 

3.6, -4.5 

3.5, -4.3 

Spiral 

0.01 

0.000 

0.000 

0.01 

0.3 

0.002 

0.3 

-0.003 

Roll 

-2.4 

-2.4 

-2.5 

-2.4 

-2.0 

-1.9 

-2.2 

-2.2 

Dutch roll 

(2.9,0.05) 

(2.7,0.08) 

(2.7,0.08) 

(2.9,0.05) 

(3.4,0.07) 

(3.0,0.13) 

(3.4,0.07) 

(3.1,0.12) 


(c) Mach - 0.9, altitude - 15,000 fL 


Phugoid (0.09,0.26) 0.0,-0.002 0.0, -0.002 (0.08,0.27) -0.005,-0.48 -0.005,-0.50 -0.005,-0.60 -0.005,-0.62 

Shortperiod 42,-65 1.7, -3.1 1.7,— 3.1 44,-69 5.1, -7.0 4.9, -6.5 49,-68 4. 7-63 

Spiral 0.01 0.0002 0.0002 0.01 0.2 0.002 0.2 -0.009 

Roll -4.9 -4.4 -4.6 -4.9 -3.4 -3.6 -4.3 -4.7 

Dutch roll (3.8,0.07) (3.7,0.10) (3.7,0.09) ( 3,8,0.07) (4,4,0.08) (3.8,0.15) (4,4,0,06) (4,0,0,14) 


(d) Mach = 1.2, altitude = 30,000 ft. 


Phugoid 
Short period 
Spiral 
Roll 

Dutch roll 


(0.07,0.22) 
2.4, -4.0 
0.006 
-4.1 

(4.2,0.08) 


0.0, -0.32 
1.7,— 2.4 
-0.0003 
-4.0 
(3.9,0.09) 


0.0,— 0.31 
1.7, -2.5 
-0 0003 
-3.9 
(3.8,0.10) 


(0.08,0.22) 
2.4, -4.1 
0.006 
-4.1 

(3.8,0.08) 


0.0, -0.40 
3.0,— 5.2 
0.208 
-3.2 

(4.6,0.08) 


0.0, -0.39 
3.0,— 4.9 
0.003 
-3.2 
(3.7,0.15) 


0.0,-037 
2.9, -5.2 
0.16 
-4.0 
(4.4,0.08) 


0.0, -0.37 
2.9,-4.9 
-0.034 
-4.0 
(3.5,0.14) 
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Table 6. Comparison of gain/phase margins between STARS models and models based 
on linearized aerodynamic data for X-29A aircraft. 


Axis 


Model 


AER07 

flexibilized 

STARS 

rigid/ 

flexibilized 

STARS 

flexible 

(a) Mach = 

0.7, altitude = 

40,000 ft. 


Pitch 




LF X/gain margin 

1.26/7.84 

1.15/6.38 

1.15/6.41 

HF X/gain margin 

23.9/11.2 

24.1/11.3 

24.2/11.4 

Gain X 

5.89 

5.53 

5.54 

Phase margin 

49.0 

48.5 

48.6 

Roll 




Frequency X 

15.2 

14.9 

15.0 

Gain X 

2.45 

2.13 

2.09 

Gain margin 

17.9 

19.2 

19.4 

Phase margin 

78.8 

74.5 

75.3 

Yaw 




Frequency X 

11.1 

11.2 

11.0 

Gain X 

3.87 

4.13 

4.06 

Gain margin 

12.0 

10.8 

10.9 

Phase margin 

62.1 

63.9 

64.6 

(b) Mach = 

0.9, altitude = 

30,000 ft. 


Pitch 




LF X/gain margin 

2.23/7.88 

2.01/6.54 

2.01/6.74 

HF X/gain margin 

23.8/6.65 

23.6/6.49 

23.5/6.47 

Gain X 

9.64 

9.16 

9.33 

Phase margin 

41.3 

41.2 

40.3 

Roll 




Frequency X 

13.7 

14.0 

14.0 

Gain X 

2.69 

2.5 

2.33 

Gain margin 

17.0 

18.9 

19.6 

Phase margin 

68.0 

57.3 

61.4 

Yaw 




Frequency X/gain margin 

13.2/11.3 

13.6/10.5 

13.6/11.2 

LG X/phase margin 

1.46/268.0 

1.47/262.0 

1.50/261.0 

HG X/phase margin 

5.64/56.1 

5.94/60.10 

5.65/62.7 



Table 6. Concluded. 


Model 

Axis 

AER07 

flexibilized 

STARS 

rigid/ 

flexibilized 

STARS 

flexible 

(c) Mach = 

0.9, altitude 

= 15,000 ft. 


Pitch 




LF X/gain margin 

2.85/7.17 

2.29/5.37 

3.20/5.76 

HF X/gain margin 

26.4/6.47 

26.1/6.22 

25.6/5.97 

Gain X 

12.3 

11.3 

11.7 

Phase margin 

38.7 

39.6 

37.5 

Roll 




Frequency X 

13.3 

12.9 

13.3 

Gain X 

3.47 

3.42 

2.93 

Gain margin 

15.0 

16.5 

18.2 

Phase margin 

66.0 

48.2 

59.9 

Yaw 




Frequency X 

13.1 

13.8 

13.7 

LG X/phase margin 

2.36/248 

2.38/241 

2.46/239 

HG X/phase margin 

5.91/54.2 

5.84/66.2 

5.45/69.4 

Gain margin 

12.6 

12.2 

13.4 

(d) Mach = 

1.2, altitude 

= 30,000 ft. 


Pitch 




LF X/gain margin 

2.08/9.98 

1.49/7.85 

2.33/8.07 

HF X/gain margin 

24.7/7.28 

24.4/7.03 

24.9/7.67 

Gain X 

9.26 

8.66 

8.17 

Phase margin 

45.8 

44.8 

43.5 

Roll 




Frequency X 

12.5 

12.5 

12.3 

Gain X 

2.62 

2.75 

2.13 

Gain margin 

16.9 

18.4 

20.4 

Phase margin 

72.0 

57.3 

72.3 

Yaw 




Frequency X 

13.2 

13.8 

13.7 

LG X/phase margin 

2.25/248 

2.21/245 

2.29/243 

HG X/phase margin 

6.73/48.2 

6.99/54.7 

6.30/60.3 

Gain margin 

10.3 

9.45 

11.0 


APPENDIX A 

FORMULATION OF THE EQUATIONS OF MOTION 

The equations of motion of an aircraft or a flight vehicle can be derived by considering various forces acting 
on the system. In classical mechanics, a Lagrangian definition of the equations enables a suitable derivation of the 
equilibrium state of a flying vehicle. This derivation requires certain physical concepts leading to the definition of 
the energy and the kinematics. The energy can be expressed in terms of kinetic energy, potential energy, and the 
virtual work done. The kinematics are generally expressed with respect to an inertial coordinate system which is 
assumed to be at rest. In the study of space-bound vehicles, such an inertial system is referenced to a distant star. 
However, in the case of an Earth-bound vehicle, the inertial system can be fixed to the center of the Earth where 

it is assumed that the angular velocities and the linear motions of the Earth have little effect on the dynamics of 
the aircraft. 


Coordinate System and Kinematics 

Let the position vector of a mass point on the aircraft be described by 

r,- = r r + r 6j + q e . (A-l) 

where 

is the position vector of the mass point i 

r r is the position vector of the rigid body-fixed coordinate system origin O with respect to inertia! space 
Tfi, is the position vector of mass point i with respect to origin O 

q ei is a vector of the deformation of the mass i with respect to the rigid body-fixed coordinate system 

If P. Q, and R are the angular velocities of the aircraft measured about the rigid body-fixed axis system, then 
the velocity of the mass can be written as 


- d, , _ _ 

r » = -q K r.) + w x r,- 


(A-2) 


where w - Pi + Qj + tffc . Equations (A-l) and (A-2) in matrix notation for all mass points can be written as 


for 


q - <t> r q r + q e (A-3) 

q = <b r q r + q e + ft(ct>rqr + q e ) (A-4) 


and 


Q 


0 

w 

0 


( A-5) 


ui 


0 -R Q 
R 0 -P 
-Q P 0 


( A-6) 


where q r and q r are rigid body displacement and velocity vectors, respectively, <t> r is a rectangular matrix of the 
rigid body mode shapes, and q e is the elastic coordinates vector. 
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Kinetic Energy 


Assuming that the mass matrix M for the aircraft system is known, the kinetic energy (KE) of the system is 
given by 


KE= 2 q M q 


( A-7) 


To simplify the algebra in the following development, the displacement vector q and the corresponding velocity 
q will be expressed in the generalized coordinate system, that is 


« = [«' O' 1 ( l ) = 


( A-8) 


where <b e is a rectangular matrix of the vibration modes. The velocity vector in terms of the generalized coordinates 


is given by 


q zr Otj + QOt] 


Substituting equation (A-9) into (A-7) gives the kinetic energy 

KE= i- + a^T)] r M [Oii+ f2<t>T)] 

Potential Energy 

If the stiffness matrix of the aircraft system is denoted by K , the potential energy of the system is 

PE=^ q r Kq = 

Virtual Work 


(A-9) 


( A-10) 


( A- 11) 


The external forces acting on the system can be classified as those dependent on the elastic deformation and the 
oscillatory motion of the aircraft, and those independent of the motion. The motion dependent forces are. 


1. rigid body aerodynamic forces, 

2. unsteady aerodynamic forces, and 


3. body forces. 

Motion independent forces are: 

1. atmospheric gust loads and 

2. engine thrust. 
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Rigid Body Air Loads 


The computed air loads in a rigid configuration involve several parameters such as angle of attack (a), angle of 
sideslip (/?), angular velocities P, Q, and R, and the control deflections 8e , 8a, and 8r related to the elevator, aileron, 
and rudder, respectively. Since the aerodynamic forces are nonlinear functions of these parameters, a set of panel 
loads are computed at discrete values of these parameters and stored in a matrix array. In the analysis the air load 
vector for derived values of the parameters is computed by interpolation. Hence the rigid air load vector (air/rigid) 
symbolically is represented by 

Fair /rigid = F (or ,/3,P,Q,R, 8e, 8a, 8r) 

The virtual work done by this force is 

VWi = t,VF ( ' air , rigid (A- 12) 

Unsteady Aerodynamic Forces 

The aerodynamic forces arising from elastic deformation and oscillatory motion can be given by 

Kir, flex = 9A e (fc)q ( A-13) 

where A e ( k) is an influence coefficient matrix and q is the dynamic pressure. The contribution to the virtual work 
from the unsteady forces is 

VW 2 = qr) T <t> T K(k)<t>V 
= <?h T [Q(*0] v 

where Q (k) are the generalized forces due to r\ r and r\ e coordinates. As k approaches 0, Q( k) corresponds to 
the steady aerodynamic loads. However the air load due to the rigid body modes such as pitch and yaw is already 
included in the F air j rigid vector. Therefore the steady air load components caused by rj r should be subtracted from 
Q ( k) to give 

0 (k) = Q (fc)- ^ -J" ( A-14) 

Then the virtual work is 

VW 2 = qri T Q(k)Ti (A-15) 


Body Forces 


Let (f>, 0, be the Euler angles defining the orientation of the body-fixed axis with respect to the Earth. Then the 
gravity vector is 

( - sin 0 'S 

cos 9 sin <{> (A-16) 

cos 6 cos <f> I 

The body forces are then given by 

Fbody - MEg (A-17) 

I 3 

where E = The virtual work due to body forces is 


VW 2 = rj"o T MEg 
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Gust Loads 


The dynamic loads in the gust environment (discrete or random) can be computed based on the unsteady aero- 
dynamic influence coefficients discussed earlier. The contribution to the virtual work from the gust loads is 

VW 4 = Tj r «I> r A e ( k)a(u) G gu3t (A-19) 

where a(w) is an equivalent downwash computed from the gust environment and G ?U3t is a measure of the 
gust intensity. 


Engine Thrust 

Future tactical fighters are expected to employ thrust vectoring to achieve superior maneuverability. There- 
fore the work done in more than one axis system is an important factor in the analysis of stability and maneuver 
performance. The virtual work contributed by the engine thrust is then 


VWs = tjVt 


( A-20) 


where T is a vector representing the components of the thrust at specified points. 


Equations of Motion 


Having defined the expressions for the energy, it is then permissible to derive the equations of motion using the 
Lagrangian principle, that is 


d_ /dL\ 
dt \dT}) 



(A-21) 


where 


L = KE - PE -VW 


( A-22) 


in which KE is the kinetic energy given by equation (A-10), PE is the potential energy given by equation (A-l 1 ), 
and VW is the virtual work as the sum of the terms given by equations (A-12), (A-15), (A-18), (A-19), and (A-20). 
Performing the indicated differentiation with respect to each component of the generalized coordinates ( 77 ) and sim- 
plifying, the following equations of motion are obtained 




gF, 


air /rigid 


ME£ + — gQcGgujt 

body force thrust 


gust 


<b r Md> 
general mass] 


Vr 

Ve 




Mfl + Q r M | <D 
\centrifugal inertia 


£K$ + g> r Q r Mfl<b +q Q(Jfc) 

general stiffness inertial stiffness unsteady 


Vr 

Ve 


Vr 

Ve 


(A-23) 


These equations have the same form as equation (4) in the main text for describing the structural dynamics, yet 
also include the additional terms to perform a nonlinear maneuver analysis. Structural damping can be included 
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with the velocity components. The inertial stiffness matrix given by O t Q t MM> is of significant importance for 
the design of the feedback control system because this term contributes to the nonlinear stiffness property of the 
control surfaces. 

The centrifugal inertia given by MO + Q^M j 3> is nonlinear in angular velocities p, q , and r , and gives 
rise to the load factor imposed on the system. The equations of motion defined by equation (A-23) can be reduced to 
a first-order differential equation and solved for the state velocities. However, for performing a maneuver analysis 
the velocity and acceleration terms (?] e and 7 )e) of the elastic modes may be assumed to be small compared to the 
corresponding rigid body terms. The second set of equations in equation (A-23) can then be solved for r\ e and 
substituted in the first set of equations resulting in the following equations of motion in the rigid body coordinates 

M % = r [qF air/rigid - ME g + F r - (M Q + Q r M ) <Mr] ( A-24) 

where r is the aeroelastic correction matrix 

r = [<t> r r + gQ re B _, O e r ] ( A-25) 

in which 

B = Kee-gQee(fc=0) ( A-26) 

The stiffness matrix K €e is an effective stiffness including the centrifugal effects. 

The expression given by equation (A-26) can be treated as an eigenvalue problem for the determination of the 
wing divergence speed (qdiv)- 

The rigid body velocity and acceleration vectors can be written as 

U 

V 

w 

Vr — p 

Q 

R 

and 

U \ 

V 

w 

p 

Q 

\ R 

Thus equation (A-24) reduces to a first-order equation in the state velocities. For a given time history of the control 
variables, a transient maneuver performance of the aircraft can be determined. 

If the inertia term M % is set to zero, equation (A-24) reduces to a nonlinear steady trim equation. 
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APPENDIX B 

INERTIAL-TO-BODY-FIXED COORDINATE TRANSFORMATION 


The inertial frame reference axes are used to compute the aerodynamic forces resulting from the free vibration 
modes of unconstrained elastic bodies. However, for stability and control analysis the preferred coordinate system 
is the body-fixed coordinate system. If the inertial frame and reference axes are inclined with respect to each other 
in a trim state through the Euler angles (¥,©,<!>), then the velocities at the aircraft center of gravity in the inertial 
frame and body-fixed reference axes are related by 


* h 


= r(¥,0,o) 



in which 


r = 


cos ¥ cos 0 cos ¥ sin 0 sin <t> — sin ¥ cos <I> cos ¥ sin © cos <t> + sin ¥ sin <t> 

sin 'F cos 0 sin *F sin 0 sin <I> + cos ¥ cos <l» sin ¥ sin 0 cos <t> — cos ¥ sin <5 

— sin 0 cos © sin O cos © cos <t> 


Similarly, the accelerations in the inertial frame are given by 

= T (¥,©,<!>) 




J B 


u 

V I + o 

w 

A x \ 

A v 

A* ) b 

where Q is the skew-symmetric angular velocity matrix. Considering the accelerations in the body coordinate system 


+ Q 


= r ('F, 0 ) «i>) 




Now with reference to a trim state at orientation ( ¥ 1 , 0j , 4>, ) , 


d> = d>i 
0 = 0 , 
¥ = 'F, 
U = Ui 
V = Vx 
W - W\ 
P = Px 

Q = Q 1 

R = R\ 
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the perturbations are 


<t> = <t>i + 

0 = ©i + e 

v = y, + 4> 

U = Ui + u 

v = y, + v 

W = W\ + w 

p = Pi + p 

Q = Qi + q 

R = R i + r 


The perturbed velocities in the inertial frame become 


forn = r('P 1 ,0 1 ,^ 1 ),and 



= r, 


Ai = 


(: 

1. 

A, | 

W \ 

x,4> 

A x,0 

A X,fj> 

y,<t> 

A y,0 

A y,'l> 

z,<t> 

A z ,9 

A Z,lj> 


The elements of Ai are as follows 


A x,* = V\ ( cos <Di sin 0! cos 'P, + sin <f>! sin 'Pi ) 

+ (cos <hi sin Vj - sin <!>! sin 0] cos T,) 

A x ,e = —U\ sin 0[ cos 'Pj + V\ sin <t>! cos 0i cos Vj 
+ W\ cos Oi cos 0] cos 'Pi 

= —Ui cos 0i sin ^i — V\ ( sin <J>i sin 0] sin *Pi + cos <t>i cos Ti) 
+ W\ ( sin <t>, cos 'Pi - cos <D] sin ©i sin ) 

A y,<t> ~ (cos <I>i sin ©i sin 'Pj — sin <J>i cos *Pi 

- W\ ( sin <t>j sin 0] sin ^ + cos O] cos qq ) 

A v,e = —U\ sin ©i sin 'P) + V\ sin dq cos 0i sin 'P, 

+ W\ cos 4>i cos ©i sin 'Pi 
Aj, ( y, = U\ cos 0i cos 'Pi 

+ Vi ( sin <J>i sin ©i cos 'P] - cos <bi sin ¥i ) 

+ W\ ( cos <hi sin 0, cos 'Pi + sin sin 'Pi ) 

A z ,<j> = V\ cos <bi cos ©i - W\ sin Oj cos ©i 
A z,8 = —U\ cos 0i - Vi sin «t>i sin 0] - W\ cos <I>i sin 0i 

A z,ij> = 0 


These equations are consistent with the formulation contained in reference 15 for general reference conditions 
with the inertia axes initially oriented to coincide with the body axes through T] . 
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The body-axis perturbation accelerations are 


Gx 

a v 

a z 


B 


u \ 

V 

w 


+ Q 



+ A 


where 


and 


Q 


A = 


0 

-R i 

Qi 

R i 

0 

-Pi 

-Q\ 

Pi 

0 

0 

W i 

-Vl 

-W 1 

0 

c/i 

Vi 

-U 1 

0 


Perturbed accelerations in the inertial frame then become 


P ^ 

9 

r / 



The state-space equations of motion in the inertial frame are given by 

x /f = A r xj r + B r u 

and the required state-space equations in the body-fixed coordinate system are obtained by the follow 
ing transformations 

X/ r = TiX Br 


(*\ 


( ^ \ 

y 


y 

z 


Z 

4> 


b 

6 


e 

i> 

= T, 

b 

X 


u 

y 


V 

z 


w 



p 

e 


9 

i> / 


r / 


where the primed coordinates denote inertial displacement components of the aircraft center of gravity projected 
onto the body coordinate system, and the relevant matrices are given as 

" Ti 0 0 O' 

0 10 0 

Tl ~ o Ai r, o 

o a 2 o r 2 
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with 


r 2 = 


COS 01 

sin <I>i 
cos 0i 


- sin q>i 
COS 01 
COS 01 


Similarly for state-space velocities 


©i tan ©i cos 
A 2 = — COS ©1 0 

elk *' ,ane ' 0 


where 


*j r = T 2 x Br + T 3 x Br 

r, ooo 
T o r 2 o o 
2 nfi n a r, o 
o o o r 2 


anu 

■ 0 Ai 0 0 

T _ 0 A 2 0 0 

3 " 0 0 0 0 

0 0 0 0 

In the following analysis no transformations are applied to elastic 
full-state vector transformations are given by 

t T i 0 ' 

1 = L 0 1 

t T 2 0 ' 

Tz = [ 0 I 

nr. _ [" T 3 0 

T3 0 0 


I 

and aerodynamic lag state vectors. Thus the 


and outputs are given by 


x/ = A xj + B u 
T 2 x B = AfiX B - T 3 x b + Bu 
T 2 x B = ^ATi — T 3 ^ x B + Bu 

= T 2 -' (Afi - T 3 )x B + T^’ftu 
= Ax b + Bu 

y = Ct,x b + Du 
= Cx B + Du 
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APPENDIX C 

STABILITY DERIVATIVES 
Determination of Static and Dynamic Stability Derivatives 

Longitudinal Derivatives 

A brief outline is given of the relations used in the evaluation of the longitudinal stability derivatives from un- 
steady generalized aerodynamic forces. For small perturbations of pitch rate, assuming a wings-level trim condition, 
(q = 6) and the rate of change of angle of attack (d), the lift and moment coefficients of an aircraft can be written 
in the Maclaurin series as 

C L = Cl„ + C La <y + Cl* (jy) + (^ v ) 

Cm = Cm 0 + Cu a oi + Cm* ( 27 ) + 

where c is the mean aerodynamic chord and V is the flight velocity. For harmonically oscillating motion at reduced 
frequency k 


O—o „«*>t 

9 ~ 6 ° e ’27 

II 

•a . 

(C-3) 

iut 

a = Q! 0 e , — 

= ika 

(C-4) 


Substituting equations (C-3) and (C-4) into (C-l) and (C-2) and omitting the constant terms Cl 0 and Cm 0 - the 
first-order lift and moment coefficients can be written as 


Cl — Cl u oi + ikC Li a + ikCi q 0 ■ ■ • (C-5) 

Cm = Cu.<* + ikCMa a + Cm • • • (C-6) 

For pure pitching motion ( 6 = a) and pure plunge (/i) 

Cl = ^f = [c La + ik(c L > + C L ,)]a (C-7) 

Cm = g gg = [ Cm* + t'fc (Cm* + Cm,)] a (C-8) 


where S is the area of the wing and Q mn is the generalized force acting on mode m due to mode n, written with 
respect to an inertial axis system aligned with the body axes (ref. 15) at a = 0. The first rigid body mode is the pure 
plunge mode and the second is the pure pitch motion. 

For a unit value of a, the stability coefficients are written as 


C La = 


Re ( OhB ) 

S 


(C-9) 
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(C-10) 


Cm, = 
Cl, + C L , = 
Cm, + Cm. = 


-Re (Qee) 

S 

ImiQhe) 

kS 

-ImiOee) 

kS 


For pure harmonic plunging motion h = h 0 e' wt , the angle of attack (a) can be written as 


(C-ll) 
(C- 12) 



9 = 0 

then from equations (C-5) and (C-6), the lift and moment coefficients for plunging motion are given by 


(C-13) 


Cl = (ikCt. - PCl.) (y) = -If ■■■ 

(C-14) 

Cm = (ikCM, - k 2 CM,) (jfj = ^ 

(C-15) 

Letting ( = 1 , the dynamic stability derivatives can be written as 


„ Re(Qhh) 

Cla ~ k 2 S 

(C-16) 

S, Re(Qoh) 

~ ~~ k 2 S 

(C-17) 

Finally from equations (C-ll) and (C-12) 


Itti(Qm) Re(Qhh) 

L '~ kS k 2 S 

(C-18) 

„ Im(Qeg) _ Re(Qgh) 

" ~ kS + k 2 S 

(C-19) 


Thus the six longitudinal stability derivatives can be determined from the pure plunge and pitching oscillatory mo- 
tions at an infinitesimally small reduced frequency. 


Lateral-Directional Derivatives 

The side force, rolling moment, and yawing moment coefficients under small perturbation assumptions can be 
represented by 

& = Cy,0 + Cy, (§) H- C,, (^) + Cy r (£)■■■ (C-20) 

a = C t ,0 + c, t (§) ♦ Q, Qfr) - ($) . (C- 21 ) 

c * = w + S (#) + c - (£) - (w) - (C - 22) 
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As in the longitudinal case, Q mn will be used to represent the generalized force acting upon mode m due to 
the motion of mode n. The rigid modes (y, <f>, i/0 are orthogonal to each other. The area of the wing is S, b is the 
span, V the velocity, and A; is a reduced frequency. The side force generated by pure yawing motion (J3 = —i>) is 
written as 


CAW = 


Rewrite equation (C-20) as 


yQvjrt 

S 


/ • ^ mQ yti i 

= -sr + '—g~* 

= s ,|Mat ( - / 5) + Sat 


<-« 


+ !& (£)<-* 


RcQyip - ImQ y $ fc 

S P kS \b 



£b_ 

2V 


(C-23) 


C y =C Y b (3 + C Yt 


0 12 V 


+ Cy. 


■(&) 

-<v+ K - »0 (If) 


Comparing equations (C-23) and (C-24), for a unit value of y one obtains 

Cy, = 
and 


Is 


0 


^ ■ (r&) (i) 


(C-24) 


(C-25) 


(C-26) 


Tb evaluate Cy. and Cy T the side force due to pure harmonic lateral motion is considered, (y = y 0 e Mt ), 


0 = V = * 


and 


then equation (C-20) becomes 


(¥) • 
p = r = 0 


(C-27) 


Cy = 


ikCy B - k 2 C Vk 


(I)] (¥) 


/2j/\ _ Qyp 

~ ~ S 


(C-28) 


Let = 1 to get 


i?e( Qyy) ( C 

Cy > = ~w 


G) 


( C-29) 
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and from equation (C-26) 


C Yt = 


ReQ 


yy 




k 2 S kS 
The side force generated by pure rolling motion (p = 0) is 


G) 


CM) = 


= -y 


= ~v 


= -y 


= -y 


S' s 

s fcS 


rfleQrt /mg rt f.-jg) J 

S' 9 kS \2vj v \ 


RzQ,y4> . ImQyf (c 

S * kS V&, 2V 



RmQyt fc 


kS 


Dm 


since ReQ v< f, is zero. Now comparing the like term in equation (C-20), and unit y, 

c ^ = -^r(D 

Hence, all four side-force derivatives are extracted from equations (C-25), (C-29), (C-30), and (C-32). 
Tlie rolling moment generated by pure yawing motion (/? = —rp) is written as 

cm) = 


= \ Re ^ + 

s s 


= <t> 

\ReQH ( _^ + ik Im ^( /?)] 

= <p 

RcQ^ IttiQ ^ 

S ( P) kS 

(&) 

2 v) 

= 4> 

ReQw ImQw fc\ ( pb\ 

S P kS \bj\2v) 


-P) 


Rewrite equation (C-21) as 


Pb 


Cl = Cl p P + Ci, [fy J + a 
= Ct p P + {c h - Cl) 


'-Pb' 
2V . 


pb_ 

2V 


(C-30) 


(C-31) 


(C-32) 


(C-33) 


(C-34) 
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Comparing equations (C-33) and (C-34), for $ one obtains 


c„ = 


and 


( ImQw 

kS 


m 


- C tr - 

From equation (C-27), equation (C-21) becomes 
Let = 1 to get 


r _ R^Qh}_ 

C e * k 2 S \b) 


and from equation (C-36) 


Ct, = 


ReQa 

k 2 S 


ImQw 

kS 


]G) 


The rolling moment generated by pure rolling motion (p = <f>) is 

CM) = 


= -4> 

s s 


= -4> 

\ ReQ »4, + ik ImQ ** J 
S kS 


= -4> 

ReQ a ^ ImQ H /.uc\ , 

. S * kS \2 v) 9 \ 

= -4> 

ReQu, , I™Qw ( # 

S * + kS {l){2V 

ImQ m (c\ ( pb\ 

= kS U)Uv) 



(C-35) 

(C-36) 

(C-37) 

(C-38) 

(C-39) 


(C-40) 


since ReQ ^ is zero. Now comparing the like terms in equations (C-21) and (C-40), and unit cf>. 

Hence, all four rolling moment derivatives are extracted from equations (C-35), (C-38), (C-39), and (C-41). 
The yawing moment generated by pure yawing motion ( f3 = —rp) is written as 


(C-41) 


c„(« - -Sfia* 


= _v» \ ReQ + i ImQ ++^\ 
s s 


= i> 

P) + /?)] 

= i> 

5 1 P kS 

(i^—) ( 

V 2VJ 

= 0 

RcQ^i/, ItyiQ $$ /c\ / 0b\ 

S P kS \b)\2v) 
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Rewrite equation (C-22) as 


Cn — Cnp/3 
= CnpP 

Comparing equations (C-42) and (C-43), for 


•«.(#) --(if) 

* («.-=.)($) 

a unit value of 0 one obtains 


c "> - \~s~) 


and 


c* - c* = 


_ / / c 


\ kS 


m 


Again from equation (C-27), equation (C-22) becomes 

a = [oo, - (?)] (|) = 

Let = 1 to get 

/nr _ ^(Q#) 

" * ~ k 2 S \b) 
ReQw ImQ± y, ] / c\ 

\b) 


and from equation (C-26) 

Cnr ~ L k 2 S kS 

The yawing moment generated by pure rolling motion (p = <j>) is 


(C-43) 

( C-44) 
(C-45) 

( C-46) 

(C-47) 

(C-48) 




= -1> 
= 

= -0 
= _Vr 


<t> + * 


5 

5 

S 

R&Q 

S 


4> + 


<t> + 

+ 


ImQ^ f.uic , , 
ImQ m (c 


kS 



<f>b 

2V 


- (?) (4) 


Its \b/\2VJ 


since ReQ^ is zero. Now comparing the like term in equation (C-22), and unit </>, 


Gv = 


IrnQ^i (c 


kS 


(S) 


(C-49) 


(C-50) 


Hence, all four yawing moment derivatives are extracted from equations (C-44), (C-47), ( C-48), and (C-50). 
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Aerodynamic Derivative Representation 


The relationship between airplane forces and moments in the inertial axis system and airplane stability deriva- 
tives is summarized in tables C-l and C-2 for the longitudinal (symmetric) and lateral-directional (antisymmetric) 
analyses, respectively. The quantity is the trim angle of attack. Column headings designate rigid body deflec- 
tion dof, and row labels represent direction of resulting air loads whether they be linear forces or angular moments 
caused by the deflections. The generalized forces at a very low reduced frequency (/c <0.01) define the generalized 
aerodynamic stiffness and damping matrices from which stability derivatives can be evaluated. If the derivatives 
are available from other sources, such as wind-tunnel or flight-test results, they may be incorporated in preference to 
the theoretically derived values. For coefficients pertaining to different reference points on the aircraft, corrections 
must be made to the derivatives (ref. 16). The inertial and body axis coordinate systems are aligned with each other 
at trim a\ . 

It is implicitly assumed in tables C-l and C-2 that the derivatives are defined with reference to the body axes. 
Longitudinal stability coefficients defined with respect to the stability axes remain unchanged, yet lateral-directional 
coefficients must be adjusted as a function of angle of attack, as in table C-3, for a proper implementation into the 
generalized aerodynamic matrices defined in the inertial axis system. 

Whether or not elastic modes are included will determine which stability derivatives (rigid or flexibilized) should 
be implemented into the matrices of tables C-l, C-2, or C-3. When all elastic modal dof (elastic generalized coordi- 
nates) are included in the dynamic analysis, the derivatives used should be only rigid coefficients (not flexibilized). 
The elastic modes included in the analysis will include the elastic increments to the rigid body derivatives through 
the generalized displacements at low values of reduced frequency. A rigid body analysis, without elastic modal dof, 
requires using the flexibilized derivatives from another source to account for elastic static deformation. No x dof 
aeroelastic forces exist since the panel aerodynamics that are not perpendicular to lifting surfaces are not computed. 
Hence, the generalized matrix coefficients for the 2 -dof generalized coordinate must include the elastic increment 
for rigid and elastic analyses. The forces corresponding to the x-dof are therefore zero for all but the lowest reduced 
frequency. When a combined rigid-elastic analysis is required with a small subset of elastic modes included, such 
as when closed-loop roots are desired in a digital system, the rigid coefficients from another source should be used. 
In all analyses, the correction rigid body data is smoothed into the real and imaginary force matrices at the lowest 
reduced frequencies. 

Flight control applications require an accurate description of stability characteristics over a frequency range 
dependent on the dynamics of the controller. In the STARS formulation, the three symmetric dof (x, z , 6) and three 
antisymmetric dof (y, <£, ^) are combined to model linear six-dof motion about a reference flight condition. Various 
methods exist to approximate flexibilized dynamics with residual approximations when some modes are deleted 
from the analysis. Presently STARS can be run with all flexible modes included or some truncated, but residual 
stiffness or flexibility effects are not used as a reduced order approximation of elastic dynamics. 
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Table C-l. Stability derivative representation in the rigid-body symmetric generalized 
aerodynamic stiffness and damping matrix elements. 

(a) Rigid body, symmetric generalized aerodynamic stiffness matrix. 


X 

Z 

e 

6 

x 0 

0 

S [Cd„ - ot\Cn„ - ai C Lmi + 

S{C Dl -*xC Ls ) 

z 0 

0 

S ( Ci a - ot\Ci v + ot\CD a — «iCd.) 

S(C Ls + axC Dt ) 

e o 

0 

Sc( Cm a + OtlC/mJ) 

ScCrn^ 


(b) Rigid body, symmetric generalized aerodynamic damping matrix. 


X 

z 

6 

6 

x (Cd. - a\ C L J 

- ai Cl 0 - Cl) 

(Cd, + Co* - ax C L , - ax Cl 

0 

z ^ (Cl. + QfiCr).) 

(Cl„ + <*i Cd„) 

JTJj" (Cl, + Cl„ + axCo, + axCoaj 

0 

0 

~&cC ma 

2^-c (-C m , - C m4 ) 

0 


Table C-2. Body axis derivative representation in the rigid-body antisymmetric generalized 
aerodynamic stiffness and damping matrix elements. 

(a) Rigid body, antisymmetric generalized aerodynamic stiffness matrix. 



V 

4> 


6 


y o 

S (-Cl — ax Cy^ q»i 

SCy, 

—SCy e 


<t> 0 

-SbCe^ax 

SbCt, 

—SbCe t 


o 

— SbCngOix 

SbCng 

-SbCm 

(b) 

Rigid body, antisymmetric generalized aerodynamic damping matrix. 


y 

4> 


V» 5 

V 

~U~x C y » 

2T7T (- C ^-aiCy A ) 

^7r(“ 

■Cy r + Cy^j 0 

4> 


^b(-C^-Ct $ a ,) 

(■ 

-a + C tfi ) 0 


-jfc bC r# 

^ b (- C ^~ C n $ ax) 

$r 6 (- 

~C Tir +C nff ) 0 


37 


Table C-3. Stability axis derivative representation in the rigid-body antisymmetric 
generalized aerodynamic stiffness and damping matrix elements. 

(a) Rigid body, antisymmetric generalized aerodynamic stiffness matrix. 



y 0 5 (-Cl - a\Cy^j c*i SCy fi ~SCy s 

<f> 0 S6 cos «i + sin ai 56 (Ce 0 cos ai — sin c*i j 56 (— Ce 6 cos ori + sin ai) 

ip 0 56 cos <*i - C^, sin c*i) 56 (Cry cos ai + Ct p sin a<i) 56 (-C^ cos a\ - Ce 6 sin aa) 


(b) Rigid body, antisymmetric generalized aerodynamic damping matrix. 


y 

<t> 

1 1 > 8 

y 

2^7T ( _Cy » “ a ‘ Cy /») 

( _Cn + Cy >) 0 


4> jjj-6 (-Cep cos on ) ^jj -6 ( — C^cos on ) 3^-6 (-C^ cos ai) 

+ ^-6(o*sin ax) sin ai ) + $T[ b ( Gnr sin cv i ) 

-^j-6 (^(cosaOa') + ^j-6 (Cg 0 cos ai) 0 




+ (c n ^(sin of] ) ofi ^ 



Cup cos 

^(-GyCOSa,) 

’fjj-b (—Cnr cos ^l) 


~~Ui^ {pfy s in <*i) 

^2T7T 6 ( C ^ sin ai ) 

~TU\b sin ai ^ 



“2^6 ( COS 

+ ^jj-6 (Cn^j COS O'! ) 0 



— ^C^(sin «i)ai) 

+ ^/rK c ^ sin ai ) 


38 



REFERENCES 


1. Gupta, K.K., STARS - A General-Purpose Finite Element Computer Program for Analysis of Engineering 
Structures , NASA RP-1129, 1984. 

2. The NASTRAN Theoretical Manual, NASA SP-22 1 (06), 1981. 

3. Noll, T., M. Blair, and J. Cerra, “ADAM-An Aeroservoelastic Analysis Method for Analog or Digital System,” 
J. Aircraft , vol. 23, no. 11, Nov. 1986, pp. 852-858. 

4. Zislin, A., E. Laurie, K. Wilkinson, and R. Goldstein, “X-29 Aeroservoelastic Analysis and Ground Test 
Validation Procedures,” AIAA-85-3091, Oct. 1985. 

5. Taylor, R.E, K.L. Miller, and R.A. Brockman, “A Procedure for Flutter Analysis of FASTOP-3 Compatible 
Mathematical Models, Vol. 1 - Theory and Application,” AFWAL-TR8 1-3063, June 1981. 

6. Gupta, K.K. and C.L. Lawson, “Development of a Block Lanczos Algorithm for Free Vibration Analysis of 
Spinning Structures,” AIAA-87-0946-CP, Apr. 1987. 

7. Abel, I., An Analytical Technique for Predicting the Characteristics of a Flexible Wing Equipped with an 
Active Flutter-Suppression System and Comparison with Wind-Tunnel Data, NASA TP- 1367, 1979. 

8. Whitbeck, R.F., D.G.J. Didaleusky, and L.G. Hofman, “Frequency Response of Digitally Controlled Systems,” 
J. Guidance Control, vol. 4, no. 4, July-Aug. 1981, pp. 423-427. 

9. Appa, K., “Constant Pressure Panel Method for Supersonic Unsteady Airload Analysis,” J. Aircraft , vol. 24, 
no. 10, Oct. 1987, pp. 696-702. 

10. Sefic, Walter J. and Cleo M. Maxwell, X-29A Technology Demonstrator Flight Test Program Overview, NASA 
TM-86809, 1986. 

11. Whitaker, A. and J. Chin, “X-29 Digital Flight Control System Design,” AGARD-CP-384, Active Control 
Systems - Review, Evaluation, and Projections, Oct. 1984. 

12. Gupta, K.K., STARS - An Integrated General Purpose Finite Element Structural, Aeroelastic, and Aeroser- 
voelastic Analysis Computer Program, NASA TM-101709, 1990. 

13. Waszak, Martin R. and David K. Schmidt, “On the Flight Dynamics of Aeroelastic Vehicles,” AIAA-86-2077, 
1986. 

14. Gera, Joseph and John T. Bosworth, Dynamic Stability and Handling Qualities Tests on a Highly Augmented, 
Statically Unstable Airplane, NASA TM-88297, 1987. 

15. Miller, R.D., R.J. Fraser, M.Y. Hirayama, and R.E. Clemmons, Equation Modifying Program L219 (EQMOD), 
Vol. 1: Engineering and Usage, NASA CR-2855, 1979. 


39 



y 



Summing junction outputs 
u * r - Gy 
X = AX + Bu 
y = CX + Du 


Figure 1. Feedback control system. 
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Figure 2. Major modules and capabilities of STARS. 
















Figure 3. STARS aeroservoelastic analysis flow chart. 
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Figure 6. Node line comparisons for G VS and STARS analyses for various modes for the X-29 A aircraft, symmetnc 
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Figure 7. Continued. 
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Figure 7. Concluded. 
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Figure 8. Unsteady aerodynamic paneling idealization (DLM and CPM). 
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Figure 10. X-29A aircraft flutter solution plots using GVS measured modal data (symmetric case). 
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Figure 11. X-29A aircraft flutter solution plots using STARS calculated modal data (antisymmetric case). 


57 




































STARS flexible 

GVS 



(a) Gain, dB. 



(b) Phase angle, deg. 


Figure 17. Longitudinal loop gains computed with GVS and STARS models for flight condition 3 
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Longitudinal loop gains computed with GVS and STARS models for flight condition 4. 
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Figure 23. Directional loop gains computed with GVS and STARS models for flight condition 1 
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Figure 24. Directional loop gains computed with GVS and STARS models for flight condition 2. 
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Flgure 27. Longitudinal closed-loop responses calculated using GVS and STARS flexible models for flight condi 
tion 1. 
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Figure 33. Longitudinal loop gains from flight data, AER07, and 
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Figure 34. Longitudinal loop gains from flight data, AER07, and STARS for flight condition 4. 
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Figure 35. Modal damping and frequency comparisons between flight measured data and STARS for symmetric 
wing first bending. 
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Figure 35. Continued. 
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Figure 35. Concluded. 
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Figure 36. Modal damping and frequency comparisons between flight measured data and STARS for fuselage 
vertical bending. 
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Figure 36. Continued. 
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Figure 36. Continued. 
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Figure 37. Modal damping and frequency comparisons between flight measured data and STARS for antisymmetric 
wing first bending. 
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Figure 37. Continued. 
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Figure 38. Modal damping and frequency comparisons between flight measured data and STARS for fuselage lateral 
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Figure 38. Continued. 
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Figure 39. Modal damping and frequency comparisons between flight measured data and STARS for fin first bending. 
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Figure 39. Continued. 
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